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Abstract 

We revisit an inequality due to Onsager, which states that the (quantum) liquid structure factor 
has an upper bound of the form (const.) x \k\, for not too large modulus of the wave vector k. This 
inequality implies the validity of the Landau criterion in the theory of superfluidity with a definite, 
nonzero critical velocity. We prove an auxiliary proposition for general Bose systems, together 
with which we arrive at a rigorous proof of the inequality for one of the very few soluble examples 
of an interacting Bose fluid, Girardeau's model. The latter proof demonstrates the importance of 
the thermodynamic limit of the structure factor, which must be taken initially at k ^ 0. It also 
substantiates very well the heuristic density functional arguments, which are also shown to hold 
exactly in the limit of large wave-lengths. We also briefly discuss which features of the proof may 
be present in higher dimensions, as well as some open problems related to superfluidity of trapped 
gases. 
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I. INTRODUCTION 



The properties of the Fourier transform of the pair (or density-density) correlation - in 
particular its extrapolation to zero angle, and the corresponding relation to the isothermal 
compressibility - have been thoroughly studied in connection with the critical behaviour of 
classical fluids (see M. E. Fisher's early review The situation is very different regarding 
the corresponding quantity - which we shall call the liquid structure factor - for quantum 
Bose fluids Q . Even for the impenetrable Bose gas - a prototype of one-dimensional systems 
with on-site repulsion - very few results on correlation functions exist, notably on the Fourier 
transform of the one-particle density matrix, a mathematical tour-de- force by Lenard jjj. 
Some results on higher-order correlations exists as well j^J, in connection with the so-called 
Fisher- Hart wig conjecture |fj, but they are restricted to Dirichlet or Neumann boundary 
conditions (b.c). 

The quantum liquid structure factor plays a crucial role in the Feynman variational 
Ansatz p, but periodic b.c. are essential in this context, because the corresponding wave- 
function must have definite momentum. 

In this paper we present a rigorous proof of Onsager's inequality for the liquid structure 
factor in the one- dimensional impenetrable Bose gas (Girardeau's model j?j|). This inequality 
implies the validity of Landau's criterion of superfluidity [8( with a definite, nonzero critical 
velocity. Some features of the proof hold - others are expected to hold - in higher dimensions, 
as discussed in the conclusion. 



II. ONSAGER'S INEQUALITY, THE LANDAU CRITERION AND THE FEYN- 
MAN ANSATZ. SUMMARY 

Several criteria for superfluidity exist in the literature [l| . By one of the standard criteria 
the free Bose gas is a superfluid, while by the Landau criterion (see, e.g. Ql), i t is not. 
Landau postulated a spectrum of elementary excitations (see, e.g., Fig. 13.11 of [10]), which 
leads, by a well-known argument, to the startling property of superfluidity of the Bose fluid, 
i.e., the fact that under certain conditions its viscosity (e.g., of liquid He-II) vanishes, i.e., 
the liquid is capable of flowing without dissipation through very thin capillaries, as long as 
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the velocity v of the liquid has an absolute value below a certain critical speed v c : 



H < v c . (i) 

In general, v c increases as the diameter of the capillary decreases. See or 11 1 as 
excellent introductory texts. 

The Hamiltonian of a system of N particles in a cube A of side L, upon which periodic 
b.c. are imposed, may be written (we use units such that h — m — 1, where m is the mass 
of the particles): 



3=1 

The momentum operator may be written 



N 



P A = -i^V r (3) 



Above, 



V = ( 4 ) 

l<i<j<N 

where $ is the interparticle potential, which we assume to satisfy the conditions appropriate 



to existence of the thermodynamic limit 



121. Let 



be the unitary operator which implements Galilean transformations. The operator imple- 
menting unitary Galilean transformation has actually the form 

However we may restrict ourselves to t — in (|5bp. and use (|5ajl instead, just for the purpose 
of obtaining the transformed Hamiltonian (energy), which indeed yields the correct formula 
(J7J). The reason for the latter is that the commutator of the generator of U# in (j5b|) with 
H\, independs of t, because 

[H A ,P A ]=0. (6) 
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The operators are defined on the Hilbert Space TL A = L 2 ym (A N ), of symmetric square inte- 
grable wave-functions over the iV-fold Cartesian product of the cube A with itself, satisfying 
periodic b.c. on A. We shall disregard the spin variable. We have 

HI = UiHjyU, = H A + v-P A + l -Nv 2 . (7) 

Suppose the fluid moves with constant velocity v with respect to some fixed inertial frame 
F . In the rest frame of the fluid the Hamiltonian is H A , and by © H A is the Hamiltonian 
with respect to the frame F. In the rest frame of the fluid, it posesses a unique ground state 
f^A, with energy E (N, L) say, and zero momentum: 



#aOy = E (N,L)Q A , (8a) 

P A n A = o. (8b) 

The fact that the ground state is simple, i.e., nondegenerate, is a property of Bosons only 
which will be used later (see, e.g., for a popular textbook account: for precise conditions 
on $ in (HJ) see where a proof is given). By (jUJ) - (JSJ), 

Hp A = (e {N, L) + ^Nv 2 J Q A . (9a) 
Thus the energy eigenvalue of Q A i n frame F is 

E* = E (N,L) + ^Nv 2 . (9b) 
Because of (jSJ, there exists common eigenstates of H A and Pa, which we denote by ip A , i.e., 



where fceA*. and 



H^l = E A (k)ip k A , (10) 

p A 4 = m (ii) 



A = { ^-n,n ee (ni,n 2 ,n 3 );ni G Z, % = 1, 2, 3 \. (12) 
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Q, 



Let ^i 0),fc be common eigenstates of H\ and Pa of the smallest energy above Eq(N, L), which 
we denote by E^k). Under certain conditions, seen to be realized in soluble models (Q, 
), it may be expected that 

E^(k) = E (N,L) + e A (k), (13) 

where e\(k) is the energy of an "elementary excitation" (see ^| for a better definition 
of this important concept, which should not be confused with "quasiparticle" - the latter 
dissipate). By (|ZJ) and JIEJ), the energy of the lowest energy excitations above the ground 
state in frame F is given by 

E (N,L) + e A (k)+v-k + ^Nv 2 . (14) 

There will be dissipation if this energy is smaller than the energy of Q\ in frame F, which 
is given by (J9b|) : thus, by (|9bj) and (JH|), there is no dissipation, i.e., there is superfluidity, if 

e A (k) + v- k> 0, (15) 

which is expected to hold if (0) is satisfied. This is Landau's criterion . It is satisfied (with 
a calculable v c ) in the two remarkable models (j?j], jlJ], I5I]). The above relation (JT5|) also 
provides the condition for Q A to be the ground state in the moving frame. 

In order to formulate Landau's criterion more precisely, we have to define the thermody- 
namic limit of e A (k) for any k, not only those in ([12)1 . Let {A n } be an increasing sequence 
of periodic boxes, 

k n G A* and k n — ► k ^ 0. (16a) 
Since any real number can be expressed as a limit of rational numbers, we may define: 

e(£)=lime An (4) (16b) 

n 

on the above defined 'e(k) we formulate then: 

Assumption i - 'e(k) > for k ^ is independent of the sequence in f!16b|) and is a 
continuous function of k when the limit on the r.h.s. of (I16bl) exists. 
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By (115(1 . (fT^|) and assumption i, the following inequality holds for all k e M 3 : 



Z(k)+v-k>0. (17) 
The above may look too pedantic, but we shall see that the thermodynamic limit will play 

— # 

an important role later on. With the above assumptions, we shall drop the tilde on e(k) 
from now on. We remark that assumption i has been proved in at least one case, that of 
magnons (or spin- waves) in quantum ferromagnets (Q|, Q|). 

An important attempt to "explain" ([T7]) . i.e., to derive (|T7(l from microscopic laws, was 
undertaken by Bogoliubov (see |8j for an excellent review), assuming the existence of Bose 
Einstein Condensation (BEC). As remarked by Leggett in his beautiful review |18], "while 
extremely suggestive, Bogoliubov's result referred to a dilute system, which is rather far 



from real-life liquid He-II" . What seems to be usually not emphasized is that Bogoliubov's 
model, even in the version which conserves particle number [8|, is not Galilean-covariant, 
i.e., does not satisfy (J7|), and, thus, does not satisfy local mass conservation, an essential 
physical requirement (this seems to be well known but, for a proof, see [19]). Bogoliubov's 
seminal work led, however, to the genesis of an understanding of e (p) (defined by (j2HI)) 
through later work (see 20( and references given there). One of its important byproducts 
was the realization that it is the repulsive part of the potential which plays the major role 
in superfluidity. This is vindicated in the remarkable one-dimensional models of Girardeau 
j]| and Lieb and Liniger (Q], ^^). The latter, however, illustrate an important fact: 
superfluidity is independent of BEC, because these models display superfluidity in the sense 
of Landau 3|> but no BEC (this was proved in [^] for Girardeau's model, and is presumably 
also true for the Lieb-Liniger model). 

n 

As remarked by Leggett in 18], a particular successful attack on the full He-II problem 



(i.e., not only the dilute case) was made by Feynman p| and Feynman-Cohen 2l( (see also 
Lieb's remarkable early review jsjj) through a variational wave-function (in our notation 

©, (HDD, CP): 



N 



^X(xi,...,x N ) = ^2e ik ^n A (x u ... } x N ) } (18) 



i=l 



(the Feynman-Cohen Ansatz j2l| leads to somewhat better results, but we shall not consider 



it in this paper). By using partial integration, and the periodic b.c. (see, e.g., ll(, Exercise 
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p. 262, or |6J), we obtain 



where 



S A (k) 



E (N,L)+£ A (k), 



(19) 



Nk 2 



2«|^ A > 25 A (fc) 



(20) 



and 



N 

1=1 



-ik-Xj 



fi|(fi, . . . ,£at), 



(21) 



(As follows from the Perron- Frobenius theorem, f2 A is unique, and may be taken to be 
positive [13]). Except for a factor p = y, (J2*T|) defines the liquid structure factor, which may 



- in case of liquid He-II - be measured independently of the excitation spectrum, by X-ray 
scattering (see Q, Fig. 13.10). 

By (|21|). <S A (&) has a direct physical interpretation: it is - for k ^ - the Fourier transform 
of the pair correlation function in the ground state ([lj], pg. 951). 

Equation (J18)) and the min-max principle imply that, for k ^ 0, (|19|) is an upper bound 
to the energies of the "elementary excitations", i.e., the lowest eigenvalues of H A which are 
larger than the ground state energy. The following assumption will be made: 

Assumption ii - The Feynman Ansatz ()18|) yields a (qualitatively) good approximation to 
the energy of the elementary excitations. 

The above assumption is widely accepted today [18| . because it is confirmed by experiment 
both for liquid He-II ^(J and for trapped gases (2^. 

Assumptions i and ii imply, together with (|20|) that 
Assumption i' - Under (jlfiaj) . the limit 



S(k) = lim<S A „(£Vi 

n 

exists and is a continuous function of k. 



(22) 
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The fact that it is necessary to take the thermodynamic limit with k ^ is shown in 
Appendix A for the free Bose gas: the limits k — > and the thermodynamic limit do 
not commute. This procedure has also a physical reason: it is necessary to eliminate the 
(infinite) forward scattering peak, i.e., by considering k 7^ 0, before probing the density- 
density correlations for the infinite system pg. 951). 

II. 1. Onsager's Inequality 

The following inequality (|2*3*j) (with condition ([24))) will be referred to as Onsager's in- 
equality Q]: 

5(jfe)<c|jfe|; MO, (23) 

where 

< c < oo (24) 

is a constant independent of N, L, but possibly dependent on p. 

We see from the above definition that Onsager's inequality, taken together with Feyn- 
man's Ansatz (Assumption ii and (fTHjl ) implies Landau's condition (JTZJ), with v c (defined 
by (1)) given by 

v c = {2c)-\ (25) 

For the free Bose gas S\(k), given by (j2U), equals 1 (see ()A2|) of Appendix A), which is, of 
course, necessary in order to be compatible with (j2*U|) . 

Let eo(p) denote the ground state energy per particle in the thermodynamic limit 

E (pL s ,L) 

e Q {p) = hm . (26) 

By jl2| e (p) is a convex function of p. Thus, when e "(p) = d 2 e (p)/dp 2 exists, 

e "(p) > 0. 



S 



We assume that 



< e "(p) < oo. (27) 
and pijl constitute a precise formulation of Onsager's inequality Q], with 

c = fey 1/2 <28 » 

Note that (J23J) and (J28J) are consistent due to (J27|) if we assume < p < oo which we do 
henceforth. 

In this paper we derive (|23jh (|24|) for a soluble model, Girardeau's model defined in 
section IV, which is believed to be a prototype of a class of one-dimensional Bose fluids 
with pointwise repulsive interactions. We believe that several aspects of the derivation are 
relevant to higher dimensions and discuss this in section V. In particular, Proposition 1 of 
section III is of general validity, as shown in Q. Our derivation in section III "rounds off" 
some points left over in Q], such as emphasizing the important role of nondegeneracy of 
the ground state - which leads to nondegenerate perturbation theory - and the vanishing of 
some mixed terms" (JSTIj) . In section IV we conclude the proof of (|23|). (}2"4"j) for Girardeau's 
model. 

It is important to realize the intrinsically nonperturbative nature (in k) of the bound (|2~H|) 
on St^(k) in the thermodynamic limit. Writing in ()21j1 . 



N 



5> 

i=i 



N 



l + ^ e -ifc-(£i-zj) ) (29) 



it becomes clear that the double sum above would contribute a term 0(N 2 )/N = 0(N) to 
S\(k) if the configurations a?, ~ x^, Vz ^ j, were not suppressed: due to the strong repulsion 
at short distances, such configurations are highly improbable in Q(x\, . . . , xn) (think of hard 
cores). But the linear term in the perturbative expansion (in k) of (|29j) is identically zero 
by symmetry, thus suggesting S\(k) = 0(|A;| 2 )! 
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III. PROOF OF AN AUXILIARY PROPOSITION 

We now provide a proof of an auxiliary proposition (Proposition 1). The derivation in 
is incomplete in some details. Let Po jA denote the projector onto the orthogonal complement 
Q A of the unique ground state Q A , V'a be given by (jTSjl . and 

A A = P 0>A (H A - E (N, L)) 1/2 P , A , (30) 
B A = P 0)A (H A - E (N, L))- 1 / 2 P , A . (31) 

bmce, by © and (HHJ), Va 

is an eigenf unction of Pa with eigenvalue k, and k ^ 0, ip A G 
Pq,\H A) and thus, by the Schwarz inequality 

<vM> 2 = (a a ^1ip a ^1) 2 < (a a ^1ia a ^)(Pa^||Pa^) (32) 

but 



(Avv||A A vi> = (p 0)A vii(^A- J Bo(^^))p OjA vi) 

= (4\(H A -E (N,L))4) 



k 2 

N- (33) 



by partial integration and use of the periodic b.c. We also have that 



(B A ^ A \B A ^ A ) = (4\(H A - E Q {N,L))- l ' 2 P^ A (H A - P (iV ) L))- 1 / 2 ^) 

= {n A \wl\H A - E Q {N,L))- l P^ A wln A ). (34) 

In flMJ), Wjlf is the multiplication operator 



N 



(W^ip)(x u ...,x N ) = ^e lk ^^{x l ,...,x N ). (35) 



i=l 



We now write = Yj%=i cos (^ ' %i) + iSi=i srn (^ • ^») hi (j^H). For finite AT, L, (P A — 
Eq(N,L)) has a nonzero lower bound on Pq,aH a , and insertion of a basis of eigenstates 
\4>i) of (i7 A — Eq(N,L)) in Po tA H A is rigorously justified because f2 A is a normalized state. 
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Furthermore, because H A is a real operator, we may choose this basis of eigenstates \<f>i) as 
consisting of real functions. The "mixed terms" in (|34|) 



N N 

-1 



^2 R _ R {^A\^sm(k ■ x^l^^jl^cosik ■ Xi)\n A } (36a) 

and 



Ei — En 

3 u i=l i=l 



N N 

^2 E _ E (fal ^2cos(k ■ Xi)\tt A )(tt A \ ^sin(fc- x^j) (36b) 

j ■? i=l i=l 

thus add to zero, and we obtain from (|34|) and analytic perturbation theory [23J: 
Proposition 1 : 

<V>M> 2 < -i-{B K ^\B^ k K ) (37a) 



2 



where 



V / A=0 



dX- 



■Eo >2 {N,L,X) 



Above, Eq } i(N, L, A) is the ground state energy of the Hamiltonian 



N 



H A 1] (A) = H A - Eq(N, L) + \J2 cos (k 



i=i 



and En t 2 (N, L, A) is the ground state energy of the Hamiltonian 



A=0 



(37b) 



(38a) 



N 



H[ 2) (A) = H A — E (N, L) + \J2 sm(fc ■ £). 



4=1 



(38b) 



By f)34|l and the argument following it the two quantities on the r.h.s. of (J37)l are also the 
second order energy terms in analytic perturbation theory (j3|, 24 1). The boson nature of 
the particles is used to the sole extent that the ground state is nondegerate ljj . The radius 
of convergence of the perturbation series in A is expected to tend to zero as (1/N), for our 
iV-particle system, but this is no source of trouble as long as A is taken to go to zero in (}3Tj) 
for fixed, finite N and L. What is remarkable in (|37|). (|38|) is the fact that the perturbation 
is an external potential. 
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IV. ONSAGER'S INEQUALITY AND GIRARDEAU'S MODEL 

A heuristic but appealing argument due to Onsager (see and appendix B) yields 

1 



lim 



L — >oo 



N (B^\\B^l) 



< ci (39) 



where 



c l = — TTTT (40) 

whenever the limit on the l.h.s. of (jSS)) exists. Putting together (J2IJ, (JS2J), (JSHJ), (jHS|) and 
P0|), we arrive at ((23), with c = (ci/2) 1 / 2 . 

Other arguments (see, e.g., (A.97)-(A.103) of [n]], "longitudinal sum rules") which yield 
S(k) ~ c\k\ for | A; | sufficiently small, all depend on the argument in Appendix B. This 
is a density-functional type of argument, which, in gene ral, is not generally justifiable for 
Fermions, whenever the ground state is degenerate [25]. For Bosons, the ground state is 
nondegenerate, and density-functional theory looks the same as for Fermions, but it has 
been rigorously established recently that, even in the dilute limit, a semiclassical description 
of Bosons is impossible (see, e.g., j^, remarks after (2.12)). Thus, the derivation of (|39jh 
()40|) given in |2| (a somewhat better version of which, with corrections, is given in appendix 
B) is conceptually open to question. But, perhaps most importantly, the derivation has a 
"miraculous" character, because in (|3~4"j) the contribution of intermediate states is expected 
to yield (const.) x due to the denominator (H\ — Eq(N, L))^ 1 and the fact that the 
lowest states are expected to have energy e\(k) ~ c\k\ above the ground state! These matters 



are dealt in a different way in the appendix of [2lj : there, the argument relies on the study 
of the terms of a perturbation theory which likely diverges. 

We now examine (I3U1) . (I3UI) on the light of one of the very few soluble models of a 

n n 

superfluid, the Girardeau model (the Lieb-Liniger model (|14j, |15J) is expected to yield 
similar results). We shall see that the results confirm ()39j) . (|4())l surprisingly well, and 
formula ([Blip , derived from density functional theory, holds exactly, in the limit k — ► 0. 



12 



We start by the Hamiltonian of the Lieb-Liniger model 



□ 



1 N 8 2 

H N,L = -oY,7^ +2d E (41) 

Z i=l ^ l<i<j<N 

with d > 0, with periodic b.c. on a box of length L. In the limit d — > oo one obtains 
a model of impenetrable bosons, with the impenetrable core shrunk to a point, previously 
treated by Girardeau Ml, and called Girardeau's model. Both models exhibit two branches 

' l ***"*iij* , *"*"^"'"7*'"'; 

seen experimentally [27| . Since the results on Girardeau's model may be obtained from the 
Lieb-Liniger model (|4*T|) by continuity as d — > oo we expect that our results are 

also true for the richer model ([41 p. and are, thus, prototypical of this class of models. 

One of the basic features of Girardeau's model is that, due to the repulsive interaction, 
the system acquires a finite compressibility (proportional to (e^p)) -1 ), in contrast to the 
free Bose gas, for which, e.g., with Dirichlet b.c, 

E (pL 3 ,L)=fpL 

where / is a constant, and thus, by pSj). e (p) = for all p. Thus eg(p) is also identically 
zero, i.e., c\ = +oo in (|3*9j) . (|4"Dj). Thus the inequality (|2*3*j) (with c given by P8|)) is trivially 
true by (jX2j). 

Performing the limit ci — > oo on ()41j) is equivalent to impose the subsidiary condition: 

if>(xi, ...,x N ) = if x j = x l , l<j<l<N. (42) 
The Bose eigenfunctions ip B satisfying (|4^|) and periodic b.c. with period L are given [t| as 



ij B = ij F A (43a) 



where 



A{xt, . . . , x N ) = Yl sgn(xj - Xi). (43b) 

3>l 

Above, sgn(x) is the algebraic sign of x, and tp F are Fermi energy eigenfunctions satisfying 
which are just the eigenfunctions of a free Fermi gas. For odd N we fill the "Fermi 
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Sphere" 



-l(N-l)<p<±(N-l) 



and the ground state energy is 

|(iV-l) 



E (N,L)= J2 ( 2 -J^) 2 = l(N-N^)(n P f. 
P =i 



Thus the ground state energy density in the thermodynamic limit equals 

1 

e (Pj = 

and thus 



e (p) = 

o 



e' '(p) = n 2 p 

which satisfies (12Tj) . 

By (EH), (E3), (SH) (-BaV'aI-Sa^a) ma y be calculated by the formula (N odd) 

i(AT-l) 

E 0>1 (N,L,\)= E ^ X ) 

P=-i(iV-l) 
i(AT-l) 

Eo t3 (N,L,\)= J2 ^> A ) 

P=-i(iV-l) 

where A) are the energy levels of the one-particle Hamiltonian 

H = H + XH X 

with 

- 1 * 
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and 



Hi = cos(kx) (53) 

on H = L 2 per [0, L], with k G A* = {^n, n E Z \ {0}} and the zero^ order level of H 
corresponding to E\(p, A) is 

E?'(P) = (54) 

with p satisfying (jUJ). The levels E 2 (p, A) are the same, just replacing Hi by H 2 = sin(fcr). 
The level p = of H is nondegenerate, all the others in (jUJ) are doubly degenerate, with 
energies (jHljl and eigenfunctions 

1 i 2?rp 1 _j 2irp N — 1 

f P ,a=i = L x , <P P ,a=2 = -^e l x ; p e {1, 2, . . . , 2 }■ (55) 



By Kato degenerate perturbation theory (|23j, 2J| (10.88)) we find the levels up to second 
order by solving the equation 

det{A(y? Pi/3 |ifi|y?p, a ) + X 2 (ip p ^\HiS p Hi\ip p ^) - {Ei(p, A) - E[ 0) )8/3 ja } = (56) 

where 

OP- 1 - 1 ? <67) 

and 

^ P (0) = 1^,1X^,11 + 1^,2X^,21. (58) 

Using (JH5J) . (|H7jl . (JHHJ) and plane- waves as intermediate states to calculate the second term 
in (jnnj), we find that equation (JHfi)) may be written 



. a — a c . 

det =0 (59) 



c a — n 



with 
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[i = Ei(p 7 A) - e{ ^ (60) 



—A 2 

-S ¥v ,k (61) 



k 2 -A{fp) 2 ^ 



2 



A A 2 k 2 

c= -S k ^ p + —5^ p>k ■ e k = —. (62) 



The discriminant of (|59|) is 4a 2 — 4(a 2 — c 2 ) = 4c 2 > 0, and thus the eigenvalues for \x are 

fi± = a ± c. (63) 

The cases p = 0, = fc, and = | are isolated values which by (|4*H|) . (j4T?j) do not affect 
the thermodynamic limit on the l.h.s. of (J3l?j) . Disregarding these, we see by (I63J) that the 
eigenvalues remain doubly degenerate and equal to a, given by (J6T| . Thus, by (jl8|) . (J49|) 
and (JEIJ), and the fact that the second term in (|4~K|) is equal to first, which may be proved, 
we find 

|(Ar-i) 

M^^M^gfllj) £ jr^Es- (64) 

P=-i(iv-D 1 L ' 

27rp /r fc 

L 2 

Clearly the r.h.s. of (|6*4*j) has a natural extension e\(k) to all A; ^ 0, obtained by interpreting 
k in (JH11) as a real variable. We now find, from ()64|) the following proposition proved in 
Appendix C: 
Proposition 2 : 

4 r p i 

e(/c) = lim ca(/c) = -Pu / r^dx. (65) 



L " 



Remark 1. The principal value on the r.h.s. of (J63j) may computed, with the result 

/TTp 2 / 27fp 

— — -rdx = — - arccoth — — ) . (66) 
_ np k 2 - Ax 2 \k\ \ \k\ 1 
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The limit k — > on the r.h.s. of ()66|) can be performed by l'Hopital's rule or more simply 
by noticing that arccoth(27rp/|fc|) « \k\/2np, so that 



f2np\ = J_ 
I i 2vr/ 



lim — — arccoth I — £ ) = . (67) 

fc^o fc| I \k\ / 27rp 



We finally have 

Theorem 1. For Girardeau's model, under assumption i' S(k) satisfies Onsager's inequality 
with 

c=^=—t (68a) 

where 



t ~ 1.0481. (68b) 

Proof. By Proposition 2, for any sequence {A n } of periodic boxes, the limit lim n eA n (/c) is 
a continuous function of k e [0, np] and it follows from appendix C that the convergence 
e\{k) —>■ e(k) is uniform in A for k in compact subsets of (0, 7rp], i.e., not containing the 
origin. Thus, by (|16a|) and assumption i', for all k 0, 



S A (k) = \imS An (k n ) < —=]im.e An (k n ) 

n V2 n 

= —=e(k) < c\k\ 
V2 



-^=lime An (Av, 
V2 n 



-^=lime An (fc) 
V2 n 



(69) 



where 



1 

c = sup 

fc£ [— Trp,TTp] I V 2 



i i 



(70) 



where t ~ 1.0481, corresponding to the value k = 7tp, by appendix D. This is Onsager's 
inequality (O, (EH), with c given bv (jfiSjl. ■ 

Remark 2. The result of density functional theory (appendix B) holds as an equality, only 
in the limit k — ► 0, because, from (ffiHj) . (JHHJ) and (jfi7j) 



lim lim 

fc— >0 iV^oo 

Zj 'OO 



N 



(npy 



(71) 
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which agrees with (|B11|) . However, from appendix D, we see that (jBll|) holds in very good 
approximation for all k. We shall comment on the reasons for that a priori unexpected 
behaviour in the conclusion. 

Remark 3. (|64|) and (J65|) show explicitly the importance of performing the limit k — > after 
the thermodynamic limit. Indeed, putting k = in (|64j) we obtain a Pv on the r.h.s. which 
diverges. More precisely, the proof of Proposition 2 given in appendix C becomes invalid. 



V. CONCLUSIONS AND OUTLOOK 

The Feynman variational ansatz is widely accepted today |18|. There is a good reason 
for that: when compared with experiment, it leads to the famous phonon-roton spectrum 
of superfluid 4 He (see, e.g., Fig. 13.11 of [10]). The ansatz is only expected to hold for not 
too large \k\ (\k\ < k c « 3.5A~ l ), becoming unstable with respect to the decay into several 
other types of excitations with lower energies It is also well established for trapped 

gases 22|. 

In this paper we assumed the validity of the ansatz (assumption iii) and attempted to 
derive a precise connection between it and Landau's criterion of superfluidity (fTTj) . This 
is achieved through the use of Onsager's inequality (|23|). ()24|). whose first part was proved 
in section III, Proposition 1, which "rounds off" the "almost rigorous proof" (see j^j) 



presented in |2|. The second part was proved in section IV, Proposition 2 and Theorem 1, 
for Girardeau's model, a prototype of a one-dimensional Bose fluid with pointwise repulsive 
interactions. Some features of this proof - in particular the importance of taking k — > 
only after the thermodynamic limit and the role of the one-site repulsion to guarantee the 
existence of the latter - will be present in higher dimensions. Indeed, as remarked by Lieb 
in the fact that " the potential is effectively a kinetic energy barrier for large 7 = d/p 
(where d is defined by (|4"T|) ) - a result that also holds in three dimensions - means that 
it is really immaterial to the particles whether they can 'ge t around each other' or merely 
'through each other'." The double spectrum predicted in 15] has been found experimentally 
(2?]], substantiating the above conjectures to some extent. 

For trapped gases, a recent remarkable proof establishes superfluidity accordin g to one of 
the standard criteria [2^. The latter is, however - as asserted by the authors of ^oj] - also 
satisfied by the free Bose gas - in contrast to the Landau criterion. The criteria are therefore 
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inequivalent and it should therefore be most interesting to prove (J2HJ), (J23j) for trapped 
gases, in the limit considered in 2(| (a review where further references can be found) by 
which the range of the potential also tends to zero. A first step is the result ( 261 ]. remarks 



after Theorem 7.1, in particular (7.2), (7.3 
n-particle reduced density matrices, and ( 



_ 7.4)) that there is 100% condensation for all 



261 ] . Corollary 7.2) (convergence of momentum 



distribution for the one-particle density matrix), but the momentum behaviour of the two- 
particle density matrix remains to be studied. 

It seems rather difficult to show (jH2), (JH|) directly for Girardeau's model: a study of 
correlation functions in the model is at present restricted to Dirichlet and Neumann b.c. 
Thus, our bound for S\(k) (with periodic b.c.) is new even for Girardeau's model. 

Finally, it is gratifying that our result (J69j) agrees exactly with ()B11|) of density functional 
theory in the limit k — > 0, because the latter is, in a sense, like the Thomas-Fermi theory, 
of semiclassical nature, which is then expected to hold in the limit of large wave-lengths. 
However, Appendix D even shows excellent agreement for all k G [— irp, irp], which may be 
intuitively understood from (|H7|) (Proposition 1) whereby the l.h.s. of (jBllj) depends only on 
the variation of certain energies with respect to the interaction with an external field. The 
latter is taken care of exactly in density functional theory (first term on the r.h.s. of (|B3D 



(El) 



and is thus a very special case, which is not affected by the problems pointed out in 
and (j3]' remarks after (2.12)). Since the arguments in Appendix B are independent of the 
dimension, this agreement is a very strong evidence that (|23|). (j24j) holds for any dimension. 
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APPENDIX A 



In this appendix we study the liquid structure factor (}2*T]) for the free Bose gas, where 



For k ^ 



By flU} 



while 



Thus, 



L 3 / 2 



L 3 / 2 



«5 A (A;) 



NL 3N J AN 
1 

1 + 



dxi ■ ■ ■ dx 



N 



N 



.\7. ;i V . — . ./ A x 



/ dxi---dx N e lk '^ 1 S]) 

l<i^j<N ' 



i | N{N-1) m 
NL 3N 

1 + (JV- 



- L^— 2 ) y ^ dx 1 dx 2 e- lk ^- S2) 



dxe 



-lkx 



S A (k) = 1 if k + 0, k e A* 



S A (0) = N if fc = 0. 



(Al) 



(A2) 



(A3) 



lim S\(0) = +oo. 



L — >oo 



(A4) 



()A2|) - (jA4|) show that it is crucial to take the thermodynamic limit with k ^ 0, and only 
thereafter the limit k — > 0, in order to find the behaviour of the liquid structure factor for 
large wave-length. Notice that ()A2jl must hold in order to be compatible with (|2T)|) in the 
free case. 
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APPENDIX B 



In this appendix we present in some detail a more precise, albeit nonrigorous, variant of 



Onsager's derivation of (J39j) . (J4Uj) . |2j. Note that there are misprints and some incorrections 
in the derivation in Q]. 

By the variational principle, for the Hamiltonian (|38aj) (((38b) is, of course, treated the 
same way), the ground state energy is 

E(N,L,X) = 'mf(if)\Hjp (\)tp) . (Bl) 

ip 

In (jBlJ) . ip are normalized test functions, for N Bosons in a periodic cube A of side L. Each 
test function corresponds to a test-one-particle density 



p(x) = N I dx 2 dx 3 ■ ■ ■ dxN\ijj(x, x 2 , • • • , %n)\ 2 - (B2) 



We compute the infimum ()B1|) in two steps 2j|: first, we fix a test-function p(x) and denote 



by {ipp} a the class of test functions with this p. Define the constrained energy minimum, 
with fixed p, as 

em = M{r p \H { i\\)r P ) 

= [ v x (x)p(x)dx + F A {p} (B3) 

J A 

where 



V\(x) = X cos(fc ■ x) (B4) 



and 



F A {p} = in£(^\H A - E {N,L)\i>«) (B5) 
is a universal functional of the density In the second stage, we find the infimum over all p: 

E (N J L,X) = mfE A {p}. (B6) 

p 

Let P\(x) correspond to the infimum in (|B6|) . and p A (x) to the infimum in (jB6|) for A = 0, 
which we assume are attained in a suitable space. 
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By the definition ()28j) of eo, it is reasonable to assume from (|B5|) . for large A, that 



Fa{p} 



eo{ P A x (x)}-e {p\x)} 



dx 



(B7) 



for some functional eo{p}, such that in the thermodynamic limit p A (x) — > p, and eo(p) is 
given by (J2ED. By (jB3j), (JEU), (B and (|B7|l : 



E {N,L,X) = / dx 



e {p£(*)}-e {p A (x)} 



+ A / dxp A (x) cos(fc • x) 



2 e o(p) y ^ 



p^(f)-p A (f) 



+ A / dxp A (x) cos(fc ■ x) 
'a 



:bs) 



for A "sufficiently large" and A "sufficiently small" . Equating the functional derivative with 
respect to p^ in (jB8J) to find the minimum, we obtain, under assumption 



4(p) 



p A (x)-p A (x) 



+ A cos(A; • x) = 



and 



p A (x)-p A (x) 



A 



— -cos(fc-f) 
eo(p) 



(B9) 



Inserting ()B9|) into ()B8|) we find for the actual minimum value 



Eq(N,L,X) = — (eg(p)) 1 / dxcos 2 (k ■ x) — A 2 (cq(p)) 1 / dxcos 2 (k ■ x) + Ap / dxcos(k-x) 



2 

A 2 



A 



(BIO) 



for fc ^ 0. Repeating the calculation for (|38bjl . we obtain the same result as the r.h.s. of 
(IB 10 J) except for the replacement of cos 2 (/c • x) by sin 2 (fc • x). Finally, for the r.h.s. of (J3HJ) 

we get 

±(B A 4\B A 4) = (e^p))- 1 ^ = -1^ (Bll) 
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from which (J3l?j) and (|4Up follow as an equality. It should be clear to the reader that the 
above derivation, while appealing is far from rigorous. Indeed, (jBll|) does not hold as an 
equality, as proved in the main text, but, as shown in Appendix D, it is very close to the 
exact result. Indeed, by (gZJ), the r.h.s. of (1HTT1) is (l/vr 2 p 2 ) = (1/tt 2 ) for p = 1. Bv (IfioT). it 
must be compared with (2/iT)f(x), where f(x) is the function plotted in appendix D. The 
value plotted for x = is (l/27r), which yields the r.h.s. of (jBll|) . but we see from that 
table that the value plotted even for x = 7r is very close to the value for x = 0. 



APPENDIX C 

Proof of Proposition 2. We must prove that the r.h.s. of ()64|) is, in the thermodynamic limit, 

8 r» dx 8 / r £ dx r p dx \ 

~ 2np V J_ np k 2 - 4x 2 ~ 27rp I J_ np k 2 - Ax 2 + J £ k 2 - Ax 2 V { } 

The difficulty lies in the fact that, in the r.h.s. of (|64|) . the "integration step" El = lixjL 
may not be directly identified with the e in $CA\ . because it depends on L. In Fig.l we show 
the graph of the function / : x — > ( kJ_ .i 
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Let us call the finite sum over the p's < on the r.h.s. of ()64j) by £j (the argument for 
the right side is the same). 
We have 



= A l + + A 2 + E 2 , (C2) 

where Y,j is the sum of the areas of all inscribed rectangles between —k/2 + and zero, 
and £ 2 is the sum of the areas of all inscribed rectangles between —up and — k/2 — e^. 
We have, by symmetry, 



while 



A x = -A 2 (C3) 



r dr 

s = / jJ^-i + «i (C4) 
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and 



dx 



7T p 



:i) 2 



X A 



(C5) 



where 



L— >oo 



for fixed k ^ 0. 



(C6) 



In order to prove (jC4|) . we note that is the sum of the areas of the regions in black in 
Fig.l which lie to the left of (— k/2 — El), while o? L is the sum of the areas in black between 
(—k/2 + El) and zero: 



f(x)dx)e L [l + o(l/L)\ 



TVfJ 



(C7) 



Now 



while 



: L j f'(x)dx 

-7Tp 



k . 
2+£L 



f(x)dx \e L [l + o(l/L)]. 



(-k/2 + x)(-k/2-x) 
1 1 



-7rp 



1 



1 



-A; - £l 



+ £l 



(-k/2-7rp)(-k/2 + 7vp) 
1 



(k/2 + irp)(irp-k/2) 



(C8) 



(C9) 



f'(x)dx 



SL 



-k/2 + x){-k/2-x) 



(V2) 2 ' -"-Jfe+Hil" (C10) 
For fc ^ 0, ((nZ| - (IHTOj) establish (jHHjl . Putting dCSJ) into Jd - JDU) we obtain that the 
thermodynamic limit of the r.h.s. of (fMJ) indeed equals /, given by (|C1J) . ■ 
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APPENDIX D 

In this appendix, we show the behaviour of the function on the r.h.s. of ()66|) : 



fix) = — arccoth 

x 




defined for x > 0. 

Setting p = 1, we construct the table: 



X 


/(*) 





0.159155 (*) 


7I-/10 


0.159288 


tt/5 


0.159689 


3tt/10 


0.160365 


2tt/5 


0.161329 


tt/2 


0.162601 


3tt/5 


0.164205 


7tt/10 


0.166178 


4tt/5 


0.168565 


9tt/10 


0.171428 


7T 


0.174850 



(*) is the limit (jfi7j). 

In Fig. 2 below, we plot f(x) for p = 1: 




0.5 1 1.5 2 2.5 3 

X 



FIG. 2: Plot of f(x) for p = 1 
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We thus see that in the interval [0, n], f is monotonically increasing, and thus the largest 
value for c, given by the r.h.s. of (|68af . is the square root of the ratio /(7r)//(0) = t 2 ~ 
0.174850/0.159155 ~ 1.0986, which equals t ~ 1.0481. 
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